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We study the equilibration and relaxation processes within the time-dependent Hartree-Fock 
approach using the Wigner distribution function. On the technical side we present a geometrically 
unrestricted framework which allows us to calculate the full six-dimensional Wigner distribution 
function. With the removal of geometrical constraints, we are now able to extend our previous phase- 
space analysis of heavy-ion collisions in the reaction plane to unrestricted mean-field simulations of 
nuclear matter on a three-dimensional Cartesian lattice. From the physical point of view we provide 
a quantitative analysis on the stopping power in TDHF. This is linked to the effect of transparency. 
For the medium-heavy 40 Ca+ 40 Ca system we examine the impact of different parametrizations of 
the Skyrme force, energy-dependence, and the significance of extra time-odd terms in the Skyrme 
functional. For the first time, transparency in TDHF is observed for a heavy system, 24 Mg+ 208 Pb. 

PACS numbers: 21.60.-n,21.60. Jz 



I. INTRODUCTION 

Time-dependent Hartree-Fock (TDHF) theory pro- 
vides a fully self-consistent mean-field approach to nu- 
clear dynamics. First employed in the late 1970's [IH1] 
the applicability of TDHF was constrained by the lim- 
ited computational power. Therefore, early applications 
treated the problem in only one spatial dimension, uti- 
lizing a very simplified parametrization of the nuclear 
interaction. Due to the increase in computational power, 
state-of-the-art TDHF calculations are now feasible in 
three-dimensional coordinate space, without any symme- 
try restrictions and using the full Skyrme interaction jSJ- 

HD]. 

In this work the Wigner distribution function [llj is 
calculated as an analysis tool to probe the phase space be- 
havior in TDHF evolution of nuclear dynamics. In com- 
parison to previous work [12) , where the Wigner analysis 
was performed in one and two dimensions, we are now 
able to carry out both the TDHF simulation and the 
phase-space analysis in three dimensions. Transforma- 
tion from coordinate-space representation to phase-space 
representation, i.e. calculating the Wigner distribution 
from the density matrix, still remains a computation- 
ally challenging problem. Here, we present a fully three- 
dimensional analysis which allows the study of relaxation 
processes simultaneously in all directions in fc-space. An 
early one-dimensional study of the Wigner function for 
TDHF can be found in [T5] . 

The paper is outlined as follows: In Sec. [H] we intro- 
duce the Wigner distribution function and discuss the nu- 
merical framework used in this work. Some benchmark 
results are presented to give an idea about the computa- 
tional resources necessary to compute the Wigner func- 
tion. We then introduce the principal observables sum- 
marizing the local or global momentum-space properties 



of the Wigner function. First the quadrupole operator in 
momentum space which gives rise to the usual deforma- 
tion parameters (3 and 7 to probe relaxation procecces 
in dynamical calculations. In addition, we define an es- 
timate for the occupied phase-space volume to obtain a 
relation between the fragment separation in momentum 
and coordinate space. 

Sec. Ill illustrates the geometrical structure of the 



Wigner distribution by means of momentum cuts for a 
few static nuclei. This is followed by a detailed discus- 
sion of the central 40 Ca+ 40 Ca collision, paying particu- 
lar attention to the effect of transparency. We discuss 
the impact of different Skyrme parametrizations on the 
relaxation behavior, as well as the dependence on the 
center-of-mass energy for a fixed Skyrme interaction. We 
also examine the influence of extra time-odd terms in the 
Skyrme functional. We complete this issue by looking 
at the 24 Mg+ 208 Pb reaction, where transparency can be 
studied for the case of a heavy system. 



II. OUTLINE OF FORMALISM 

A. Solution of the TDHF equations 

The TDHF equations are solved on a three- 
dimensional Cartesian lattice with a typical mesh spac- 
ing of 1 fm. The initial setup of a dynamic calculation 
needs a static Hartree-Fock run, whereby the station- 
ary ground states of the two fragments are computed 
with the damped-gradient iteration algorithm |14[ I15j . 
The TDHF runs are initialized with energies above the 
Coulomb barrier at some large but finite separation. The 
two ions are boosted with velocities obtained by assum- 
ing that the two nuclei arrive at this initial separation 
on a Coulomb trajectory. The time propagation is man- 
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aged by utilizing a Taylor-series expansion of the time- 
evolution operator [16] up to sixth order with a time step 
oft = 0.2 fm/c. The spatial derivatives are calculated us- 
ing the fast Fourier transforms (FFT). 



B. Computing the Wigner function 



The Wigner distribution function is obtained by a par- 



tial Fourier transform of the density matrix p(r- 
with respect to the relative coordinate s = r 



,t), 



/w( r .M) "/ (27r) 3 



V(r--,r+-,t), (1) 



(2) 



Because /\y is not positive definite, it is misleading to 
consider the Wigner function as a phase-space probability 
distribution. We will refer to the appearance of negative 
values for /w in Sec. |III| 

Evaluating the Wigner function in six-dimensional 
phase space is still a computational challenge and only 
possible employing full Open MP parallelization and ex- 
tensive use of FFT's. The determing factor is the grid 
size, which results in 



Nl log (Ay * N 2 y log (Ay * Ng log (N z 



(3) 



steps to provide the Wigner transform in full space, where 
N x , N y , N z are the grid points in each direction. 

Figure [I] shows the time taken to evaluate the Wigner 
function for one single time-step on a 24 2 x 36 grid for a 
16 0+ 16 collision. Benchmarks are shown for two differ- 
ent CPU's. Storing the Wigner function reduced to the 
reaction plane, i.e. ■ (x,y = 0, z, k) at one time step 
will consume ~ 140 Mb of disk space for the presented 
case in Fig. [T] Going to larger grid sizes, needed for heav- 
ier systems, and/or storing the full three-dimensional 
Wigner function will clearly result in entering the Gb 
regime. 



C. Observables 

In this section we discuss some of the observables used 
in our analysis. In order to avoid any misunderstand- 
ings we will label all observables evaluated in momentum 
space with a subscript k, and all observables in coordi- 
nate space with a subscript r. 



1. 



Quadrupole in momentum space 



As an observable to probe relaxation in phase-space 
quantitatively, we evaluate the quadrupole operator in 
momentum space. The local deviation of the momentum 
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FIG. 1. (color online) Wall clock time t in seconds (a) and 
speed-up (b) by evaluating the Wigner function at one single 
time step on a 24 2 x 36 grid for 16 0+ 16 as a function of the 
number of used processors n. CPU's used: (CPU): Intel(R) 
Core(TM)2 Quad CPU 2.66GHz, Cores: 4 (red), Intel(R) 
Xeon(R) CPU X5680 3.33GHz, Cores: 12 (blue). 



distribution from a spherical shape is a direct measure 
for equilibration. The local quadrupole tensor in /c-space 
is given by 

Q%(r,t) = y'd 3 fc[3(fc l (r,t))(fc J (r,t))-(k 2 (r,t))^] , 

. (4) 

using the m-th moment from the local momentum distri- 
bution 

with (k(r, t)) denoting the average local flow 
Jd 3 fck/^(r,k,t) 



(k(r,i)) 



(6) 



/d3fc/W(r,k,t) 

The spherical quadrupole moments Q|°(r,i) and 
Qf 2 (r,t) are computed by diagonalization of Q % £(r,t) 



Qnr,t) = \/— A 3 



96vr 



(A 2 - Ax) 



(7) 
(8) 



with A3 > A 2 > Ai labeling the eigenvalues of Q^(r,t). 
Switching to polar notation the observables 

(3 k (r, t) = v //3| (r,t) + 2/3 2 2 2 (r,t) (9) 

V2p 22 (r,t) 180°, 
7fc r, t = arctan — — — — , 10 
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FIG. 2. (color online) Two-dimensional k x -k z subspace from the full Wigner distribution function /' 3 ^(r,k) in the spatial 
center of the static nuclei 16 (a), 24 Mg (b), 40 Ca (c), and 208 Pb (d). One-dimensional cuts along the k x /k z -ax\s are projected 
at the walls. The Skyrme interaction SLy6 was used for all cases. 



are obtained via the dimensionless quantities 

,20/ ^ 47tQ 2 



/?f(r,i) 



5r^(r,i) 

4ttQ 2 2 
5r?/>(r,i) 



where 



r k (r,t) = y/(k(r,t))y P (r,t) 



(11) 
(12) 

(13) 



accounts for the local rms-radius in fc-space. The norm 
is defined such that 



p(r,t)= J d 3 fc/^(r,k,t). 



(14) 



In the presented formalism it is straightforward to de- 
fine global observables. The global quadrupole tensor is 
calculated by spatial integration 



#(*) = J d 3 rp(r,t)Q^(r,t). (15) 



Applying the same diagonalization as in the local case i7] 
we end up with a global definition for f3u°(t) and /3^ 2 (i 



For the following results we will mainly use the global 
definition since it is more compact and allows the simul- 
taneous visualization of multiple time-dependent observ- 
ables. We will however show some local results in Sec. 

nun 



2. Quadrupole in coordinate space 

To illustrate the global development of a reaction, we 
will also use the expectation value = (Qr°) of the 
quadrupole operator in coordinate space. 



3. Occupied phase space volume 

To give a rough measure for the phase-space volume 
occupied by the fragments during a heavy-ion collision 
we assume a spherical shape of the local momentum dis- 
tribution. Adding up the fc-spheres 



V k (r,t) 



4-7T 



;k 2 M)> 



3/2 



(16) 



160 MeV 



leads to the total occupied phase-space volume 



V k {t)= / d 3 r14(r,t). (17) 



III. RESULTS AND DISCUSSION 

It is the aim of this work to provide a quantitative 
analysis of the magnitude of relaxation processes occur- 
ring in TDHF. Therefore we will vary a single reaction 
parameter, while all the other parameters are fixed. The 
40 Ca+ 40 Ca-system provides a suitable test case. The oc- 
currence of transparency will also be examined for the 
heavy 24 Mg+ 208 Pb system. However, we first start with 
a brief discussion concerning the geometrical structure of 
the full Wigner distribution for static nuclei. 



A. Static nuclei 

In this section we present slices through the full six- 
dimensional Wigner distribution for some static nuclei. 
For this purpose the Wigner distribution is plotted at 
the spatial center, and for fixed momentum in the y- 
direction, i.e. k y = 0. 

Starting with 16 O, which corresponds to subplot (a) 
in Fig. |2j the Wigner function reveals strong negative 
values at the center. The hole in the momentum dis- 
tribution, which indicates the influence of distinct shell 
effects, was already visible in the lower-dimensional anal- 
ysis [12], however in that case it corresponded to small 
positive values. 

Considering heavier nuclei, we observe the appearance 
of two holes for case (b), reflecting the deformation of 
24 Mg in coordinate space. In contrast we find a pro- 
nounced peak for 40 Ca, case (c). The momentum distri- 
bution for 208 Pb in subplot (d) reveals again a central 
hole which this time is positive. This agrees with the 
reduction in the spatial density observed for 208 Pb. 

Although the geometrical structure differs from nu- 
cleus to nucleus, it is important to note that for cold 
ground-state nuclei the momentum distribution differs 
considerably from a Fermi distribution, even for a heavy 
nucleus. 



B. 



3 Ca+ 40 Ca 



We now investigate the relaxation process in momen- 
tum space for heavy-ion reactions. We have chosen the 
40 Ca+ 40 Ca-system as a benchmark to investigate the 
quantitative effect of some parameters on the relaxation 
process. All calculations in this section were done for cen- 
tral collisions (impact parameter b — 0). The numerical 
grid was set up with 36 x 24 2 grid points. 
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FIG. 3. Global observables (3 k (t) (a), y k (t) (b), V k (t) (c), 
and Qf{t) (d) are shown for a central 40 Ca+ 40 Ca collision 
with a center-of-mass energy E c , m . = 160 MeV. Each curve 
corresponds to a different Skyrme force as indicated in the 
legend. 
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TABLE I. Effective masses of Skyrme parametrizations used 
in this work are listed in connection with the maximal /3k- 
values from the plots (a, e) in Fig. [3] 



Skyrme force 


rn * /m 




SkM 


0.79 


0.0116 


SLy4 


0.70 


0.0111 


SLy6 


0.69 


0.0106 


SkI4 


0.65 


0.0102 


SkI3 


0.57 


0.0095 



1. Variation of the Skyrme force 

In the first set of calculations we vary the Skyrme 
parametrization. Figure [3] shows the results of a central 
40 Ca+ 40 Ca collision with the Skyrme parametrizations 
SLy4, SLy6 pi], SkMs [IS], SkI3, and SkM [19]. While 
SkMs was chosen as an example for an outdated interac- 
tion, the SLy(X) set of forces was originally developed to 
study isotopic trends in neutron rich nuclei and neutron 
matter with applications in astrophysics. The Skl(X) 
forces take the freedom of an isovector spin-orbit force 
into account. This results in an improved description of 
isotopic shifts of r.m.s. radii in neutron-rich Pb isotopes. 

The global development of the reaction is visualized 
in subplot (d). The time-dependent expectation value 
Qr°(t) shows the five trajectories initially in good agree- 
ment but finally fanning out. A similar splitting behavior- 
depending on the employed Skyrme parametrization was 
already found in [20] . While the two Sk(X)-forces show 
a full separation of the two fragments, there is a slight 
remaining contact between the fragments for the case of 
SLy6, which will result in complete separation in a longer 
calculation. However the trajectories for SLy4 and SkMs 
show a merged system in the final state, which was found 
to persist in long-time simulations. 

We now consider the relationship between the observed 
characteristics in coordinate space with the dynamics in 
phase space. Subplot (a) shows the /?fc-value, measur- 
ing the global deviation of the momentum distribution 
from a sphere. The initial /3fc-peak is strongly damped 
for all five Skyrme-forces. While the time development 
for all parametrizations remain in phase up to the second 
peak, later it starts to vary and continue with damped 
oscillations. For a better visualization the first peak is 
magnified in subplot (e). The taller the /3fc-peak the 
longer the fragments will stick together in coordinate 
space. The effect appears to depend on the effective mass 
m * lm. Smaller effective masses give rise to a smaller 
/^-deformation. Table [I] summarizes the m * / m- values 
associated with the maximal deformations (3™ ax for all 
the Skyrme forces used in this work. Subplot (b) shows 
the 7^ -value which indicates, whether a deformation is 
prolate, oblate, or triaxial |21j . For the present scenario 
of a central collision the 7fc-value jumps between prolate 
and oblate configurations indicating that the momentum 
distribution oscillates between being aligned primarily in 




z[fm] 

Video 1| (color online) Two-dimensional z-fc z -slice from 
the full six-dimensional Wigner distribution fff(r,lt) for a 
central 40 Ca+ 40 Ca collision with a center-of-mass energy of 
£c. m . = 160 MeV. 



the beam direction or transverse to it. For the sake of 
completeness we additionally present the occupied phase- 
space volume (c) which will prove more useful for the next 
reaction parameter to be discussed: the center-of-mass 
energy. 



2. Variation with the center-of-mass energy 

As a second reaction parameter the center-of-mass en- 
ergy, -Ec.m., is varied. Results are presented for ener- 
gies ranging from E c . m . = 2 McV/nucleon up to E c , m . — 
3 MeV/nucleon. The Skyrme interaction now is fixed to 
be SkM. For the case of the lowest (highest) energy 
Video [T] (Video [2]) provides a video visualizing the re- 
action in phase space. The calculation done with the 
lowest energy -B c . m . = 160 MeV shows two fully sepa- 
rated fragments in the exit channel. In contrast, the case 
with the highest energy (as well as the one at an inter- 
mediate energy) results in a merged system. The global 
observables are plotted in Fig. [4] It may not be obvious 
at first why the fragments should split for lower ener- 
gies and merge for higher ones. But the estimate for the 
occupied phase-space volume Vk{t) presented in subplot 
(c) indicates that V increases with energy. Therefore the 
fragments' average distance in phase space is larger, while 
in compensation they can come closer to each other in co- 
ordinate space. However, this behavior is also dependent 
on the particular Skyrme force used and the presence of 
time-odd terms, which is discussed in the next subsec- 
tion. 
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FIG. 4. Global observables /3*(t) (a), j k (t) (b), V k (t) (c), 
and Qf{t) (d) are shown for a central 40 Ca+ 40 Ca collision 
with fixed Skyrme interaction SkI4. Each curve corresponds 
to a different center-of-mass energy. 



t = O.Ofm/c 
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Video 2 (color online) Same as Video [2] with a center-of- 
mass energy of E c . m . = 240 MeV. 



3. Influence of time- odd terms 

Skyrme energy-density functionals are calibrated to 
ground state properties of even-even nuclei [TTrHTj] . This 
leaves the choice of the time-odd terms in the functional 
(current j 2 , spin-density s 2 , spin kinetic energy density 
T, and the spin-current pseudotensor J) largely unspec- 
ified [23]. Galilean invariance requires at least some of 
these terms to be present depending on the presence of 
the associated time-even term, e.g. j 2 for the (pr — j 2 ) 
combination. In our calculations we always include the 
time-odd part of the spin-orbit interaction. In order to in- 
vestigate the effects of the remaining time-odd terms, we 
have compared different choices by using a single Skyrme 
parametrization and the same test case. We choose the 
force SLy4 and start with the minimum number of time- 
odd terms which is needed to ensure Galilean invari- 
ance [23]. In the next stage, we include also the spin- 
density terms proportional to s 2 . Finally, we also add 
the combination which includes the tensor spin-current 
term (s ■ T — J 2 ). As shown in Fig. [5] at least for the 
quantities /3fe and Q 20 , varying these time-odd terms has 
a very little effect in the initial contact phase and the 
dynamical behavior becomes somewhat different only in 
later stages of the collision. On the other hand, small dif- 
ferences near the threshold energy (the highest collision 
energy for a head-on collision that results in a composite 
system. At higher energies the nuclei go through each 
other) can have large long-term effects on the outcome 
of the collision. For example a small difference in dissi- 
pation may be enough to influence the decision between 
re-seperation or forming a composite system. We have 
also checked a broader range of collision energies from 
the fusion regime up to deep inelastic collisions. The in- 
teresting quantity is the loss of fragment kinetic energy 
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between the entrance and exit channels. It was found 
that the spin terms contribute small changes to this loss 
which can go in both directions, less dissipation near fu- 
sion threshold and more dissipation above. Subplot (b) of 
Fig. [5] shows the effects near the Coulomb barrier where 
spin terms reduce dissipation, 
(a) 
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FIG. 5. Global observables /3 k (t) (a), and Qf\t) (d) are 
shown for a central 40 Ca+ 40 Ca collision with fixed Skyrme 
interaction SLy4. 
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FIG. 6. (color online) Mass number A and charge Z of 
the cluster emitted in the forward direction (feed-through) 
for central 24 Mg-|- 208 Pb reactions as a function of the center- 
of-mass energy Ec.m.. 
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Video 3, (color online) Two-dimensional z-fc z -slice from 



the full six-dimensional Wigner distribution /^'(r,k) for a 
central 24 Mg+ 208 Pb collision with a center-of-mass energy of 



E c . n 



350 MeV. 



C. 24 Mg+ 208 Pb 

The study of the 24 Mg+ 208 Pb reaction is motivated by 
experiments indicating some evaporation of heavy nu- 
cleon clusters around zero degree relative to the beam 
direction in the 25 Mg+ 206 Pb collision [22]. The follow- 
ing results are obtained on a 32 2 x 48 grid with Skyrme 
interaction SLy6. 

In the following we define as "feed-through" the mass 
A and charge Z numbers of the nuclear matter cluster 
emitted in the forward direction. Fig. [6] shows TDHF 
results on the feed-through for central 24 Mg+ 208 Pb reac- 
tions. However, the density of the emerging fragments is 



very low. We identify the energy threshold for a possible 
feed-through at E c . m . w 350 MeV. 

Video [3] provides a video showing the central 
24 Mg+ 208 Pb reaction slightly above the feed-through 
threshold in the Wigner picture. As time elapses 24 Mg is 
absorbed into the phase-space volume of 208 Pb. We ob- 
serve a rotation of phase-space density inside the merged 
system. After a short time a jet is visible, leaving the 
24 Mg+ 208 Pb compound. This corresponds to low-density 
fragment in coordinate space with a total mass compa- 
rable to an a— particle for the actual energy. 

In Fig. [7] we compare the global observables for a cen- 
tral and a peripheral reaction with impact parameter 
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FIG. 7. Global observables /3 k {t) (a), j k (i) (b), V k (t) (c), 
and Qr°(t) (c) are shown for a 24 Mg+ 2()8 Pb collision. The full 
curve corresponds to a central reaction and the dashed curve 
to a non-central reaction with impact parameter b = 4 fm. 
Both cases are calculated with Skyrme force SLy6 and with 
a center-of-mass energy E c , m , = 350 MeV. 



b — 4 fm. Subplot (a) shows that the /3/c-peaks are di- 
minished for the case of the non-central collision. In this 
case we find that the alignment of the momentum dis- 
tribution probed via the 7fc-angle (b) no longer shows a 
sharp transition between a prolate and an oblate config- 
uration. The smoother development of 7^ is due to the 
possibility for triaxial configurations in the non-central 
collision. 
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Video 4 (color online) The local observable fik (r , t) is plotted 
in the reaction plane, i.e. r = (x, y — 0, z) for a 24 Mg+ 208 Pb 
collision. The calculation is done with Skyrme force SLy6 and 
a center-of-mass energy of E c . m . = 350 MeV. 



The time development of the local observable /3&( r ) 
can be studied in Video [4] At t — fm/c we observe a 
non vanishing local j3 k deformation in the 24 Mg fragment. 
As we have previously shown in Fig. [2] (b), 24 Mg reveals 
a ground state deformation in coordinate as well as in 
momentum-space. In the reaction plane a domain of high 
j3k values occurs while 24 Mg is penetrating through the 
208 Pb fragment. The initial disturbance travels through 
the lead fragment while it exhibits a strong decrease in 
magnitude. When the low density nuclear cluster leaves 
the merged system, the /3k excitation has completely sub- 
sided. 



IV. SUMMARY 

We have presented a geometrically unrestricted frame- 
work to study nuclear dynamics within TDHF in the full 
six-dimensional phase space. The impact of different re- 
action parameters on the outcome of a heavy-ion collision 
was studied in detail for 40 Ca+ 40 Ca and 24 Mg+ 208 Pb 
systems. We find that there is transparency in both 
collisions, which is clearly reflected in the global asym- 
metry of the Wigner momentum distribution. The sur- 
prising result that in some cases the system merges at 
higher energies and shows transparency at lower ones 
can be related to the interplay between momentum- and 
configuration-space volumes which is a reflection of the 
Pauli principle. It is also interesting that the two dis- 
tributions in phase-space never truly combine to form a 
single distribution. This clearly indicates that two-body 
collisions will be necessary to achieve true equilibrium as 
the reaction proceeds to longer contact times. The de- 
tailed degree of relaxation found depends on energy and 
also the properties of the Skyrme force, where especially 
the effective mass seems to be important. The presence 
of additional time-odd terms in the Skyrme functional 
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appears to have a complex impact on the outcome of a 
collision as well. In this paper only one non-central col- 
lision was studied. A systematic investigation of impact 
parameter and energy dependence as well as even heavier 
systems would be highly interesting but is beyond com- 
putational feasibility at the moment. 
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Appendix A: The full Skyrme functional 

Following the convention used in Ref. |25j we can write 
down the full Skyrme energy density functional £skyrme- 
It depends on seven local densities and currents, namely 
the spatial density p q (r) (time-even), the kinetic density 
Tq(r) (time-even), the current density j g (r) (time-odd), 
the spin density s q (r) (time-odd), the spin-current den- 
sity </ 0>M! ,(r) (time-even), and the tensor kinetic density 
F g>Jlt (r) (time-odd). They are defined as 



s <?( r 

F*A* 
with 



P,(r,r')| r=r , 

s g (r,r')| r=r , 

V- V'p,(r,r')| r=r , 

V- V s 9 >,r')| r=r , 

-|(V-V')p g (r,r')| r=r , 
-^-Vj.K^r.r')!^ 



(Al) 



cr 

s g (r,r')= ^(nr, r'a'){a'\a\a). (A2) 

crcr' 

We will consider the recoupled forms of the proton and 
neutron densities to isoscalar and isovector densities. The 
recoupling of the density pq(v), for example, yields 



Po( T ) = Pn( T ) + PpO) > 

= Pn(r) ~ p p (r) . 



The Skyrme full functional now reads 

Skyrme = / d 3 r ]T {<#[#>] p\ + C s t \p ] S 
J t=0,l 

+C t % t A /9t + Cf s (V-s t ) 2 
+Cf s s t -As t + C;( Pt T t -j 2 t ) 

Z 

1 z 



+C, F 



fl — X 



fi.f — X 



+C t v ' J (p t V- J f + s f •Vxj t ) 



(A4) 



The functional accounts for the central, tensor, and spin- 
orbit interaction. All possible bilinear terms up to second 
order in the derivatives are included. The coupling con- 
stants of the central and tensor part are listed below in 
terms of the well known Skyrme parameters 



(A3) 



cg = 


a 4s 




C{ = 


1,1 , 1 A 
4 v 2 1 24 v 2 


+ ^3 


Q = 


1 A , 1 A 
-4*0(2 -^o)- 34*3(3 


- ^3 


C{ = 






cs = 


3 1 ,5 s 
— U + -t 2 (- + x 2 ) 
16 4 V 4 ; 




c\ = 


1/1 \ 1/1 

-g*i(2+*i) + ^(2- 




cl = 


1/1 \ 1/1 


+ a: 2 ) 


c? = 


1 11/ 

U H to (t e 

16 16 8 y 


-to) 


o — 


9 1/5, 
^64* 1+ 16* 2 (4 + ^ 




c* p = 


3 A s 1 A 
32^2 +X1 > + 32^2 


+ £2) 


^<As 


3 A ,1/1 
32*^2 ~ Xl > + 32 h k 




cf s = 


— ii H 1 2 (3i e + 

64 1 64 2 32 v e 




Bf = 


-^(t e + 3< ) 






3, 

-g(*e-t ) 




nVs 

Bq - 






bY s = 


~(3*e + *o). 





3 

32 1 



(A5) 
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The coupling constants of the spin-orbit part of the func- 
tional can be rewritten into the form 

CY J = -\w . (A6) 
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